Assisted inflation in Bianchi VIo cosmologies 

J. M. Aguirregabiria 1 '", P. Labraga 1,6 and Ruth Lazkoz 2 ' c 

1 Fisika Teorikoaren eta Zientziaren Historiaren Saila 

Euskal Herriko Unibertsitatea, 644 Posta Kutxatila 

48080 Bilbao, Spain 

pg \ 2 Fisika eta Matematika Aplikatuaren Saila 

Deustuko Unibertsitatea, 1 Posta Kutxatila 
48080 Bilbao, Spam 

^_ . a wtpagagj @lg . ehu.es, wtblalop@lg.ehu.es, c rlazkoz@eside. deusto.es 

February 7, 2008 
> 

o 
o 

r — ■ Exact models for Bianchi VIo spacetimes with multiple scalar fields 

C^) with exponential potentials have been derived and analysed. It has been 

shown that these solutions, when they exist, attract neighbouring solu- 
tions in the two cases corresponding to interacting and non-interacting 
f\ ' fields. Unlike the results obtained in a previous work dealing with the late- 

PjH. time inflationary behaviour of Bianchi VIo cosmologies, the knowledge of 

I ' exact solutions has made possible to study in detail the occurrence of in- 

Kfy nation before the asymptotic regime. As happened in preceding works, 

here as well inflation is more likely to happen with a higher number of 
non-interacting fields or a lower number of interacting scalar fields. 



Abstract 



d • 1 Introduction 



The cosmological inflationary scenario, as pioneerely proposed by Guth M, has 
been much favoured due to its troubleshooting capabilities 0. Furthermore, it 
has recently gained the support of observational evidence jjj| Q]. It is usually 
assumed that an epoch of accelerated expansion in the Universe was driven 
by one scalar field slowly rolling down its potential. Even though several scalar 
fields should enter the general picture, it is commonplace to make the additional 
assumption that only one of them played a significant dynamical role. 

Recently, Liddle, Mazumdar and Schunck p] put forward an alternative 
model representing a slight departure from those assumptions. By means of 
exact FRW examples, they showed that multiple scalar fields with exponential 
potentials can assist each other in the realisation of inflation, even if the indi- 
vidual fields are not flat enough to produce inflation on their own. This is the 
reason why they called this behaviour assisted inflation. 



The work by Liddle and collaborators was generalised by Copeland, Mazum- 
dar and Nunes along two different lines. On the one hand, they introduced 
cross-coupling terms in the potential and showed that this seems to hinder infla- 
tion rather than assist it. On the other hand, they analysed cases with multiple 
exponential potentials and showed that, by choosing the slopes adequately, the 
expansion rate can be augmented. Motivation for the use of exponential po- 
tentials is can be found in dimensionally reduced supergravity theories (see for 
instance ||, J7J, and the references therein). 

Assisted inflation presents remarkable novel features from the dynamical 
point of view. One would normally expect that a single scalar field should 
eventually dominate the dynamics, but peculiarly inflationary solutions turn 
out to be late-time attractors. This fact has been reflected in several studies 
concerned with the dynamics of assisted inflation ||, g, [To] . 

Kanti and Olive |0, [12[ explored another interesting aspect of the problem 
using a realisation of the assisted inflation proposal based on the compactifi- 
cation of a fifth-dimensional Kaluza-Klein theory. In particular, they showed 
that assisted inflation can be a remedy to the initial conditions problem in the 
chaotic inflation scenario |l3|. Along this line, Liddle and Kaloper |L4j have 
found that in this kind of inflation the spacetime does retain some memory of 
the conditions that gave rise to it, even though a homogeneous and spatially 
flat universe is generated. 

Although most of the references in the literature are concerned with FRW 
models, there is no reason to assume that at the onset of inflation the universe 
was as isotropic as it seems to be today. In fact, it is commonly believed that it 
was precisely an epoch of accelerated expansion that isotropised the universe. 
This is the reason why it is of much relevance to extend the studies on assisted 
inflation to Bianchi spacetimes. 

Our starting point is a recent collaboration of one of us |l5f , which focused 
on the study of assisted inflation at late times in FRW and Bianchi I cosmologies 
by using exact solutions and in Bianchi VI spacetimes by means of asymptotic 
approximate solutions. We construct here a family of exact Bianchi VIo space- 
times with exponential potentials which include the exact counterparts of the 
mentioned approximate solutions. By resorting to these exact, yet particular, 
solutions we have been able to perform a detailed analysis of the occurrence of 
inflation in this framework. 

Let us recall that if the models behave as desirable they will undergo accel- 
erated expansion in their early epochs, preferably for a limited period of time. 
Clearly, asymptotic studies, as those usually done up to the date, have the 
drawback of being unable to determine whether that occurs. 

As will be shown in the following sections, we have found the necessary and 
sufficient conditions to be fulfilled by the slopes of the individual potentials for 
the expansion to become accelerated. In the light of the results, we have been 
able to illustrate the cumulative effect of non-interacting fields in the occurrence 
of inflation to which assisted inflation owes its name. What is more, our results 
exemplify as well the hindering effect of direct interactions between the fields. 
For the sake of completeness we have also included an analysis showing the 



asymptotic stability of the solutions as well as an study of the isotropization 
of our models. Summarising, our results confirms previous studies and give 
additional support to the assisted inflation proposal, providing particular exact 
solutions that attract neighbouring solutions and show inflation at the start of 
their life span. 

2 Exact solutions to the multiple scalar field 
problem in a Bianchi VIq spacetime 

In what follows we are going to consider the problem of n scalar fields minimally 
coupled to gravity [J and driven by exponential potentials in a Bianchi VIo space- 
time. First, it will be assumed that the scalar fields interact through a product 
of exponential potentials; then, we will turn to consider the alternative case of 
uncoupled fields evolving via individual exponential potentials. Nevertheless, it 
should be pointed out that, although we will be loosely speaking of interacting 
and non interacting models, the geometry will be responsible for some inter- 
action among the fields even in the examples devoid of direct couplings in the 
potential. 

Before we enter the details of the two different cases to be studied, let us 
recall the form of the line element of a Bianchi VIo spacetime: 

ds 2 = _ e f(t)( dt 2 _ dz 2) + G ^ e z dx 2 + e '-dy 2 ). (1) 

For convenience we define two vectors belonging to an n-dimensional Euclidean 
space as follows: 

4>= (01,02, ■ ■ ■,</>«), (2) 

k= (fci,fc 2 ,...,/c„). (3) 

Moreover, we also demand k be constant with respect to an orthonormal basis 
of the space. 

2.1 Non-interacting fields 

In these models we will take the potential responsible for the interaction among 
the fields to be of the form 

n 

V($)=J2V oi e- k '+*. (4) 

Motivated by the fact that all scalar fields in these kind of configurations seem to 
tend to a common limit H , we will make the simplifying assumption of having 

1 Even though most of the papers on the subject deal with models in which the fields 
are minimally coupled to gravity, inflation can also be assisted in the case of non-minimally 
coupled multiple scalar fields with exponential potentials |l(|. 



4>i = (f>2 = • • • = 4>n = <t>, fei = ki = . . . = k n = k and V i = V o2 = ■ ■ ■ = V on = 
V . Additional motivation for this assumption is that it simplifies considerably 
the task of finding exact solutions. 

It will be possible to determine the model's features upon the resolution of 
the system formed by the Einstein equations 

"' = Wo- 2 (5) 



and the Klein- Gordon equation 



p--n°' v =°- m 



A first integral to (J7]) is given by 



kG m 



2nG nG' 



(8) 



where m is an integration constant. It is straightforward to realise that by 
means of this first integral and eqs. (||) and (g) the problem reduces to finding 
the solution to a single ODE, namely 

/I 1.2 \ 1 2 

GG 2 - GGG +[---—) GG 2 + -GG 2 + — G = 0. (9) 

\2 An) 2 n 

Only a particular solution to the case m — of eq. (0) has been known to the 
date J17j . Remarkably, it also admits a rather simple family of solutions that has 
not been noticed before. In order to simplify the expression of these solutions 
we introduce a new constant defined as 



Vk 2 ~2n 
m , (10) 

Vk 2 + 2n 



Up to time-shifts and rescalings of the spatial coordinates x and y, the mentioned 
solution to eq. (g) is given by 



„ fJ , A 2 \ ■ , V2n(t + to) /l 2 \ , V2n(t + t ) ,„, 

G{t) = (-2 +p ) smh ^w=wr + {-2- p ) cosh ^7w^' (11) 

and it exists provided that k 2 > 2n. Note that the solutions with p ^ have as 
their late time limit the p = case, the already known solution. 

At this point, the complete determination of the remaining metric function 
/ and the scalar field <j> is a straightforward task. Nevertheless, we will omit the 
corresponding explicit expressions because they are rather lengthy and of little 
use for the forthcoming discussion. 



The freedom to choose the value to will allows us to avoid unwanted features 
in the solutions, in particular we will choose it carefully so that we prevent 
signature changes in the p = cases. To this end we set: 

Vk 2 -2n 2p 2 -l 

t = =— arctanh (12) 

\J2n 2p z + l 

and thus guarantee G > and G > in the whole range we are going to consider 
from now on: t > 0. Note that no signature changes arise in the p — case, so 
there is in principle no restrictions upon the choice of time origin and we will 
set to = just out of utter convenience. A consequence of these choices is that 
unless p y£ the curvature scalars will blow up at the beginning of times, i.e. 
there will be a spacelike singularity of big-bang type. 

2.2 Interacting fields 

In these other models the potential responsible for the interaction among the 
fields will be of the form 

n 

V($) = l[V ol e- k ^. (13) 

j=i 

The system will then be completely determined upon the resolution of the set 
formed by the Einstein equations 



/ G 

e j — 



2VG' 



(14) 



^-i(^V--/+- = -« (15) 

G 2\G G J 2 v v v ' 



and the Klein- Gordon equation 



G '-> i 

-(/)-e f Vk = 0. (16) 

G 



One can easily check that the vector 

r G r m 

^2G k+ G (1?) 

is a first integral to (urn and rn another constant vector belonging to the n- 
dimensional Euclidean space, that is, a set of integration constants. Again, by 
using this first integral and eqs. (|1J) and ( |l5|) the problem reduces to finding 
the solution to a single ODE, which does not differ much from the one for the 
previous case: 

/i 1.2 \ .. . 1 .. 

GG 2 - GGG + f - - — J GG 2 + -GG 2 + m 2 G = 0. (18) 



It follows that the solutions to eq. ( fijl ) can be obtained from those to eq. (g) by 
remembering that now k 2 = \k\ 2 and m 2 = \rh\ 2 and performing the replacement 
n — ► 1. Correspondingly, we introduce a new constant vector defined as 



Vk 2 



V2 + fc2 



We have then 
G{t) = 



P 



sinh 



V2(t + t ) 

VF - 2 



1 

2 

IpI 



p 



cosh 



V2(t + to) 



(19) 



(20) 



Vfc2^2 
Clearly, such solutions exist 



where we have used the shorthand p 2 
provided that k 2 > 2. 

The arguments concerning possible signature changes given for the models 
in the previous section apply to these other ones too, so we set 



to = 



Vk 2 



V2 



^arctanh-^ 



1 



2p 2 



1 



(21) 



if p ^ 0, and io = otherwise. 

Summarising, the two families of spacetimes we have constructed are quite 
analogous in their expressions. It remains to see, however, if these similarities 
can be extended to their behaviour. 



3 Asymptotic behaviour and stability 

The asymptotic behaviour of the kind of ODE leading to our two families of 
Bianchi VIo geometries is worth studying, as this analysis will enable us to 
discuss their stability. In terms of the variable f2 defined by 



» = §■ 



eqs. (||) and ( |l8| ) become 



n 



n + K 



i 

2h? 



~G 2 



(Q + l)h = 0. 



Here, we have introduced a couple of new parameters defined as 

K 

and 



*1 

An 



M z = — 
n 



(22) 
(23) 

(24) 

(25) 
(26) 



in the non-interacting case. The corresponding expressions for the interacting 
case are obtained by the replacement rules given above. 

Equation (p3) has the fixed point solution f2i = — 1. In addition, if G — » oo 
asymptotically, then 0,2 — is also a fixed point solution. They correspond to 
G ex t and G oc e t '^ 2K respectively. Now, in order to determine whether these 
two solutions are stable on the G — ► oo regime we expand each of these solutions 
about fixed points by making either il = Oi + e or O = Oa + e with |e| -C 1. For 
Qi we get 

«-se* (27 > 

and it indicates that Oi is unstable. In fact, the corresponding solution G oc £ 
is spurious because it does not satisfy the Einstein equations (cf. eqs. (04) and 
(pp). The case ^2 is more interesting because it corresponds precisely to the 
late-time approximation of the exact solutions given above, in this case we get 

e=-^=e, (28) 

and this proves the asymptotic stability of our solutions for K > 0. Therefore, 
for a given m the solution f^ is the attractor for all the solutions that are close 
to it. 

4 Inflation and other kinematic aspects 

It is well known that those spacetimes for which the gradient of the scalar 
field is timelike can be reinterpreted as perfect fluid induced geometries. By 
taking advantage of this customary reinterpretation we will be able to study the 
kinematic features of our spacetimes. For our purposes we just need to know the 
expression of the four velocity of the fluid in terms of the scalar field's gradient, 
namely 

u a = ^4=nf a,/3 = 0,...,3. (29) 

In order to find out whether the model inflates or not it is necessary to look at 
the sign of the deceleration parameter q = —8~ 2 (3 9 }a u a + 9 2 ) where the scalar 
9 = u a . a is the expansion of the fluid |lqj . In the case of interacting fields, and 
after some algebra, we have arrived at the following expression, valid for any 
spacetime obtained from (|20|): 

oy + c 2 G + c 3 G 2 
9 = 2— -2—, (30) 

( c 4 + c 5 GJ 

where 

ci = Ap 2 [k 2 (k 2 + 2n)+6n 2 (k 2 + 4 n)] , (31) 



c 2 = 4 \p\ \k\ sign(p k) (k 2 + 4n) Vfc 4 -4n 2 , (32) 

c 3 = (k 2 -2n) 2 (k 2 +4n), (33) 

c 4 = 2 |p| |fc| sign(p fc) y/k 2 + 2n, (34) 

c 5 = yjk 2 -2n (k 2 + 4n). (35) 

The equivalent expression for the interacting case could be obtained by making 
the replacements \p\ — ► p = |jf|, |fc| — ► fe = |fc|, n — ► 1 and sign(pfc) — > cos A, 
where 

fc • to fc • p . . 

cosA=- = - r ^-, 36) 

fcTO Kp 

and using 

d = 4p 2 [k 2 (fc 2 + 2) cos 2 A + 6 (k 2 + 4)] (37) 

instead of (pl|). Safely, the denominator of g never vanishes for i > 0, as can be 
easily tested. 

Our first step in this analysis is to investigate the late time behaviour of q. 
We get 

2fe 2 -4n 

hm q = —5 — , (38) 

which is positive in the range of validity of our solutions, thus indicating that 
our models do not inflate at late times. This is a very useful result indeed and 
will be used later. 

Let us try to find out now whether there is a period of time for which q 
becomes negative, recalling that this is the condition for the existence of a 
period of accelerated inflation. A quick look at the expression for q shows that 
sign(pfc) and cos A < are respectively the necessary conditions for inflation, 
since otherwise q would be definite positive. In general, the quantity q will either 
vanish at a couple of time instants which we will denote t— and t + , or will never 
become null at all. Moreover, from the fact that at late times q is necessarily 
positive it follows that q will be negative between those hypothetical t_ and t + . 
The question to formulate is whether those instants exist and whether the sign 
change of the deceleration factor at t+, which ends inflation, occurs during the 
model's history t > 0. 

For our purposes it suffices to give the implicit equations of i_ and t + , in 
the interacting case these read 

G(t±) 

(fc 2 -2n) 5 v Vfc'+4H 

' (39) 



(40) 



2 bl ( V 6 V ki n - k2 ( fc2 - 2 ) n2 - 


-2fc 2 n 3 + 8n 4 


(P-2nP Vk 2 +4n 




\k sign(pfc)vfc 2 + 2n), 
solved provided that 




k 2 l9n-l 

- + 1<\/ T- 

n V n — 1 





Typically, in multifield models and for a given k, a minimum number of fields 
will be needed for inflation to occur in multifield models. Note that in the 
particular case n = 1 there will be inflation for any value of k. Let us assume 
now that we are dealing with a multifield case and that condition (E3) is fulfilled. 

Clearly, inflation will occur if t+ > and, under the assumptions made 
before, we can conclude that there will be an epoch of accelerated inflation if 
G(t + ) > G(0). Remarkably, this condition follows from (Eo|), so the latter is 
actually not only a necessary but also a sufficient condition for the existence 
of inflation. With these results at hand we conclude that our non-interacting 
models are examples of the assisted inflation phenomenon. As a matter of fact, 
in the cases in which the individual potentials are too steep inflation can still 
occur thanks to the cooperation of the fields if there are enough of them. 

The implicit equations of t— and t+ in the interacting case are slightly dif- 
ferent: 

2d / V&\/8-k 2 (P + 2)sin 2 A , , 

G(t±) = *-— r ± i ■ ky/k 2 + 2 cos A (41) 

V ' (fc 2 -2)i I v/F+4 

and they can be solved provided that 



l<\- r r. (42) 

V 1 — cos z A 

By recalling that in this case k 2 = k 2 + . . . + k 2 we can see that the more fields 
there are, the more difficult it becomes for ( fJ2| ) to be satisfied and the less likely 
inflation is. Note that in the particular case cos A = — 1 there will be inflation 
regardless of the value of k. In general lines, the role played by the number of 
fields in also agreement with existing results. The standard interpreted of this 
behaviour is a direct consequence of the way the friction depends on the field 
population. Moreover, if ( f42[ ) holds, then G(t + ) > G(0) is fulfilled too, and we 
draw the conclusion that the fulfillment ([42]) is in fact a necessary and sufficient 
condition for inflation. 

Finally, it is interesting to look at the late time degree of isotropy of the 
spacetimes under discussion. According to an intuitive criterion for a perfect 
fluid model to become isotropic at late times |19f the quotient between the shear 
and the expansion of its fluid must tend to zero on that very limit. We use here 
the standard definition for the shear scalar of the fluid [H , namely 

<y' 2 = \° a pa afi (43) 

where a a p is the shear tensor. In general, in the non-interacting case one has 
for our geometries 

a fG-G 



V3(/G + 2G) 
and at late times it becomes 



(44) 



a k 2 -2n 

hm - = — = . (45) 

*^oo6» V3(fc 2 +4n) V ; 



The latter limit is a monotonically increasing function of k 2 /n. So the larger the 
number of non-interacting fields, the more isotropy at late times. Conversely, 
using n = 1 and k 2 = k\ + ■ ■ ■ + k\ we see that a larger number of interacting 
fields will result asymptotically in a more anisotropic spacetime. 

5 Conclusions 

Families of exact solutions to the set of Einstein-Klein-Gordon equations for 
Bianchi VIo geometries with multiple scalar fields with exponential potentials 
have been derived. The models obtained include cases with both interaction and 
non-interaction among the fields, and we have shown the asymptotic stability 
of them all. In particular this means that those spacetimes characterised by 
parameters compatible with assisted inflation act as attractors of other possible 
solutions with the same values of those parameters. 

Performing a kinematic analysis based on exact, yet particular, examples 
like the ones used here has a clear advantage in comparison with studies just 
concerned with asymptotic features. The severe limitation of not having any 
information about early stages in the evolution of the model is overcome here, 
and therefore we can formulate and answer the more natural question of whether 
inflation exists soon after the beginning of times. 

A thorough analysis of the conditions under which inflation occurs has 
yielded the result that in the interacting cases profusion in the number of fields 
redounds to the improbability of inflation; whereas in the cases without interac- 
tion all the contrary happens, this is, inflation is assisted. A discussion on the 
late time isotropisation of these spacetimes has been included as well. 

In our opinion, it would be of interest to extend this work by performing an 
analogous in-depth analysis about the conditions for inflation in the setup of 
other isotropic and anisotropic cosmologies. 
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